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Abstract. We study equilibrium configurations of non-Euclidean plates, in which 
the reference metric is uniaxially periodic. This work is motivated by recent 
experiments on thin sheets of composite thermally responsive gels [1]. Such sheets 
bend perpendicularly to the periodic axis in order to alleviate the metric discrepancy. 
For abruptly varying metrics, we identify multiple scaling regimes with different 
power law dependences of the elastic energy £ and the axial curvature K on the sheet's 
thickness h. In the h — > limit the equilibrium configuration tends to an isometric 
embedding of the reference metric, and £ ~ h 2 . Two intermediate asymptotic regimes 
emerge in between the buckling threshold and the h — > limit, in which the energy 
scales either like h A/5 or like h 2/3 . We believe that this system exemplifies a much more 
general phenomenon, in which the thickness of the sheet induces a cutoff length scale 
below which finer structures of the metric cannot be observed. When the reference 
metric consists of several separated length scales, a decrease of the sheet's thickness 
results in a sequence of conformational changes, as finer properties of the reference 
metric are revealed. 
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1. Introduction 

Thin elastic sheets are being used extensively in both natural and manufactured 
structures to form elaborate three-dimensional (3D) configurations. Complex 
configurations are often a result of the material's internal structure, rather than due 
to the application of external forces or constraints. In organic tissues, for example, 
pattern formation is often a consequence of the tissue's nonuniform growth, or a 
response to differential swelling. 
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The connection between growth, or swelling, and the resultant 3D configurations 
is highly non-trivial. Different formulations of effective elastic theories of growing 
sheets have been proposed [2, 3, 4]. In most approaches, gradients of growth (or 
swelling) across the thickness of the sheet induce a non-zero reference curvature 
tensor, whereas lateral gradients of growth endow the sheet with a 2D (generally 
non-Euclidean) reference metric [5, 6, 2, 7, 8, 9]. 

Every configuration adopted by the sheet has associated metric and curvature 
tensors (also known as first and second fundamental forms [10]). Deviations of the 
actual metric from the reference metric cost stretching energy, whereas deviations of 
the actual curvature from the reference curvature cost bending energy. Shape selection 
is determined by an interplay between stretching and bending energies; the selected 
configuration (i.e. the actual curvature and metric tensors) is the one that minimizes 
the total elastic energy. 

Recent developments of experimental techniques [7, 11, 1] have promoted the 
joint theoretical and experimental study of pattern formation in numerous systems. In 
particular, various biologically-inspired elastic problems have been studied recently 
using formulations of the type described above. Metric-driven shape selection has 
been studied in the context of leaves [12], fungi [13], petals [14], and various ribbon- 
like configurations [8, 15]. Curvature-driven shaping was shown to dominate shape 
transitions in flowers [16], pines [17], and seed pod opening [18]. A similar geometric 
approach gave rise also to reduced theories for non-Euclidean rod-like structures [19], 
and to a quantitative study of shape transitions in awns [20] and isolated cells [21, 22]. 

In all the above theoretical studies, the reference metric and curvature tensors 
were assumed to be smoothly varying relative to the lateral dimensions of the sheet. In 
reality, biological materials are often highly heterogeneous, both in their elastic moduli 
and in their growth/swelling profiles. The smooth reference tensorial fields that were 
used in the models have therefore to be viewed as homogenizations of fluctuating fields. 
Homogenization was successfully applied, for example, in [23], were a "digital" 
metric, of shrinking dots embedded within a less metrically responsive matrix was 
replaced by a smooth effective reference metric. Similarly, the fibrous structure of 
a pod valve [18] and a cell wall [22] were successfully replaced by smooth effective 
reference metrics and curvatures. 

Naively, one would believe that the homogenization of the reference tensors 
is applicable when the thickness of the sheet is large compared to the scale of 
metric fluctuations. Yet, the precise interplay between small scale metric fluctuations 
and the thickness of the sheet remains to be investigated, both experimentally and 
theoretically. In particular, it is not clear what to expect when the reference metric 
consists of several, well separated, length scales. Such questions are of major 
importance for the modeling of biological tissues and for the modeling of responsive 
composite materials, and are the main motivation for the present work. 

This paper studies a common biological structure — a monolayer of a fibrous 
tissue. Fibrous tissues are ubiquitous, and appear in various sizes, ranging from 
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cell walls of plants [24] and bacteria, and up to macrostructures. Fibrous tissues 
typically consist of aligned fibers embedded within a matrix. In some cases, such 
as in sclerenchymal tissue (rigid supporting plant tissue) [25], the swelling/growth 
properties of the matrix and the fibers are different. As a result, changes in 
environmental conditions, such as humidity, or the active growth of the organ, lead 
to the buildup of internal stresses that are followed by shape transitions. 

It was recently shown [26, 18] that sheets made of two super-imposed fibrous 
layers behave like shells (i.e., have a non-zero reference curvature), and undergo 
various shape transitions. More recent experiments studied the response of monolayers 
of fibrous tissues to environmental changes. From a theoretical point of view, such 
structures form a new kind of non-Euclidean plate, with a reference metric that has 
a striped structure; it is homogeneous along the fibers, and inhomogeneous albeit 
periodic in the perpendicular direction. Such non-Euclidean plates have not been 
studied previously, and one would tend at first to think that such surfaces are flat 
"on average". It was therefore of a surprise when these monolayers were found to 
bend along the fibers. It is one of the goals of the present paper to resolve this pattern 
formation mechanism. 

The first question of interest is why would a monolayer that is for all practical 
purposes homogeneous across the layer bend? We show that bending is the natural 
response to a uniaxially periodic non-Euclidean metric, i.e., it is a response to a metric 
incompatibility. 

Another characteristic of the experimental system studied in [1] is an almost- 
discontinuous reference metric. We show that as a result it is no longer possible to 
partition the range of parameters into "thick sheets" and "thin sheets"; one has to 
account for an interplay between the thickness of the sheet, the scale over which the 
reference metric varies between bulk values, and the macroscopic dimensions of the 
stripes. In particular, we identify up to five distinct regimes of thickness, in each of 
which the the dependence of the configuration on the parameters is different. One 
of these regimes was not previously known, and its discovery sheds light onto the 
experimental and numerical results reported in [11]. 

These results show that the homogenization of heterogenous materials is indeed 
highly non-trivial and scale dependent. Moreover, we expect that the multiplicity of 
energetic and conformational regimes could be even richer in systems in which the 
reference metric consists of several well-separated length scales. 

2. The model 

We model the composite sheets studied in [1] using the formalism of incompatible 
elasticity [27, 6], and more specifically the dimensionally-reduced model of non- 
Euclidean plates [2]. A non-Euclidean plate is a two-dimensional surface, endowed 
with a Riemannian reference metric a. Other parameters of the plate are an elastic 
modulus Y, a Poisson ratio v, and a thickness h; all three parameters may be, in 
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general, position dependent, however in this work we assume that they are spatially 
constant (note, however, that in the system studied in [1] the two gels have very 
distinct elastic moduli). 

We start by setting the reference metric a. For striped sheets, it is convenient to 
choose a system of coordinates with one of the axes parallel to the stripes. We denote 
points on the sheet by coordinates (u, v), with the constant-z; parametric lines parallel 
to the stripes (see Figure IB) The symmetry of the system implies that the reference 
metric depends only on the u-coordinate. Moreover, it is always possible to choose a 
so-called isothermal, or conformal parametrization, such that the reference metric has 
the form 

' 1 \ 



a{u,v)=f\v)\ 



1 



(1) 



where f(v) is a local swelling profile (shrinkage, if / < 1). 
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Figure 1. Illustration of the experimental setting and observations in [1]. (A) Gel 
stripes that swell at different ratios upon stimulation arranged in an alternating array. 
(B) A surface composed of alternating stripes. Upon stimulation internal stresses 
cause the respective compression and tension of the two types of gels (illustrated by 
yellow arrows). (C) Functional dependence of the swelling profile on the v coordinate. 
The three curves represent swelling profiles with different transition widths, A. (D) 
Cartoon of isometric embedding in the case of an abrupt metric transition. The sheet 
curves along the stripes while undulating in the perpendicular direction. The lines 
that are at the top of the undulations are longer than lines that are at the bottom of the 
undulations. 

For simplicity, we assume that prior to the stimulation both gels have the same 
width w, so that f(v) is periodic with period 2iv . Inspired by the experimental setting, 
we assume that / ranges between bulk values 1 ± A/2, with a transition layer of width 
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A«a;. Specifically, we choose the following swelling profile (Figure 1C), 

A taring" 1 sin(f )) 

W= X+ 2 tanhA- < 2 > 
but all our results are insensitive to the precise functional form of f(v). 

Our goal is to study the equilibrium configuration of an unconstrained plate 
endowed with the metric (1). We recall that any surface is uniquely determined 
(modulo rigid transformations) by its two fundamental forms: the first fundamental 
form, or metric a and the second fundamental form, of curvature tensor b. In a non- 
Euclidean plate, the elastic energy associated with a configuration is the sum of 
two contributions: a stretching energy, which penalizes for deviations of the metric a 
from the reference metric a, and a bending energy, which penalizes for deviations of 
the curvature tensor b from a reference curvature tensor b = (non-zero reference 
curvature tensors distinguish plates from shells [19]). 

Naively, one would think that a state of zero energy can be attained by adopting a 
configuration in which a = a and b = b. This is however not always possible, because 
the fundamental forms a and b are not independent; they are constrained by a set of 
differential relations — the Gauss-Mainardi-Codazzi equations [10]. Non-Euclidean 
plates and shells differ from Euclidean ones in that the reference forms a and b do not 
satisfy these constraints, and as a result, the equilibrium configuration is selected by 
a competitive interplay between both stretching and bending energies. 

The precise form of the stretching and bending energies are material dependent. 
For small metric strains and small curvatures, it is possible to approximate these 
energetic contributions by terms that are quadratic in both strain and curvature. 
Specifically, the total elastic energy £ may be approximated by 



£ 



(W s + W B ) ds, (3) 
where Q is the domain of parametrization, Ws is the stretching energy density, 

W s = ^A a ^(a-a) a ^a-a) y5 , (4) 
and Wb is the bending energy density, 

W B = ^A a ^ 8 b aj} b y6 (5) 

(as standard we adopt Einstein's summation convention). Here A is an elastic modulus 
tensor, which for locally isotropic materials takes the form 

A «m = _L_ + a^a^) . (6) 

l+v\l-v / 

We take the domain of integration Q to be 
Q = [0J] x [-w,w], 

with periodic boundary conditions in v; i is the length of the stripes, which we assume 
to be longer than any other length scale, and ds denotes the differential area element, 

ds = ^Idetadudv. (7) 
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To summarize, the equilibrium configuration of the plate is postulated to be the 
(^-periodic) minimizer of the total elastic energy (3), where W$ and W B are given by 
(4) and (5), and the reference metric a is given by (1) with swelling profile / given by 
(2); the important parameters of / are the swelling factor A, and the transition width A. 

3. Analysis 

3.1. Why do thin striped plates buckle? 

Experiments show that thin striped plates of the type introduced in the previous 
section buckle, i.e., they adopt a bent equilibrium configuration that breaks the plate's 
planar symmetry. In this section we explain in qualitative terms the mechanism that 
leads to this buckling. A more thorough analysis follows in the next subsections. 

The pref actors in front of the stretching and bending energy densities (4),(5) have 
different dependences on the thickness of the plate h: the prefactor of the stretching 
content is linear in h whereas the prefactor of bending content is cubic in h. For 
thick plates, the bending energy is dominant, and therefore the energy minimizing 
configuration is planar. For thin plates, the stretching energy is dominant, and the 
plate will buckle, or bend, if it can thus lower its stretching energy. 

In a flat configuration the striped plates are metrically strained, because a flat 
configuration implies that both types of stripes have the same length, whereas 
their reference lengths are different (Figure 1A,B). Bending parallel to the stripes' 
longitudinal axis, thus turning the sheet into a cylindrical envelope, does not change 
this metric mismatch. If however we superimpose on top of a longitudinal bending 
a transversal undulation that is correlated with the stripes, then stripes that are at the 
top of the undulation are longer than stripes that are at the bottom (Figure ID) and 
therefore the strain is reduced. 

Such deformations cost bending energy, so that the equilibrium configuration has 
to be selected by an interplay between the stretching term, which as just described 
favors bent configurations, and the bending term, which favors planar ones. In [1], 
a simplified semi-quantitative model was derived to predict shape selection in the 
extreme case of stiff fibers embedded in a soft swelling matrix. In the next section 
we perform a quantitative analysis that reveals a fine balance between the stretching 
and the bending energies. We show that as the thickness of the sheet is reduced, the 
equilibrium configuration may be characterized by up to five distinct regimes, each 
characterized by a different dependence of the system's parameters. 

3.2. Surfaces of revolution 

We turn to study the minimization of the elastic energy (3), with a swelling profile 
(2). While a full solution of this problem is still beyond reach, a much simplified 
problem is obtained if restricting the space of configurations to surfaces of revolution — 
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an ansatz that is natural given that surfaces of revolution are consistently observed 
in experiments. 

For surfaces of revolution, both first and second fundamental forms are diagonal 
(in the chosen parametrization) and only depend on the transverse coordinate v; we 
adopt the standard notation 



«= n n\ b = n x, ■ (») 





The four periodic functions E(v), G(v), L(v), and N(v) are constrained by the Gauss- 
Mainardi-Codazzi equations [10], which given the imposed symmetries reduce to a 
pair of differential-algebraic equations, 

, E' E' E' 1 E'G E" 

L = — L + —N EN=— + — -. 9 

2E 2G 4E 4G 2 w 

By straightforward integration methods we may obtain L 2 in terms E and G, 

L 2 = cE-g, (10) 

where c is an integration constant. A direct substitution yields N 2 as well in terms of 
E and G. 

We have thus reduced the problem to that of minimizing the elastic energy £ 
expressed as a functional of two 2w-periodic functions, E(v) and G(v), and a single 
constant of integration c. Since the problem is symmetric in the longitudinal direction 
we may minimize the energy per unit length, which amounts to taking t = 1; also we 
set for concreteness the Poisson ratio to be zero. Under these choices, the total energy 
measured in units of the stretching modulus Yh/8 takes the form, 

£ = £ ((£ - ff + (G - ff) * 4 £ (L 2 + N*) d y 2 . (ID 

The minimization of the energy (11) with respect to E(v), G(v), and c can be performed 
easily by numerical optimization methods. 



3.3. Numerical resuts 

Computational results are shown in Figure 2-4. In Figure 2 we plot equilibrium 
configurations for a choice of parameters of A = 0.1, A = 1/160, and three different 
values of h. These surfaces reveal two main features: (i) the smaller h is, the larger is 
the longitudinal curvature; (ii) the smaller h is, the more noticeable is the small scale 
of the swelling profile. Both features reflect the approach of the actual metric a to the 
reference metric a. 

Figure 3 A displays a log-log plot of the elastic energy (11) as function of the 
thickness h f or A = 0.1, and different values of the transition width A (the values are 
specified in the legend). 

Four distinct regimes are observed: 
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Figure 2. Equilibrium configurations of one periodic unit for A = 0.1, A = 1/160, and 
three difference values of h. The smaller h, the larger is the longitudinal curvature, 
and the more noticeable is the metric transition. 



(i) Planar configuration: For thick sheets we obtain, as expected, an energy that 
does not depend on h, which characterizes the plane-stress solution. The buckling 
threshold h% is almost independent of the transition width A. 

(ii) Near threshold: The energy dependence on h near the buckling threshold is 
almost independent of the transition width A. 

(iii) Isometric configuration: For h — > 0, we have £ ~ h 2 , which indicates that the 
equilibrium configuration converges to an exact isometry of the 2D reference 
metric a. Note however the strong dependence of the prefactor of h 2 on the 
transition width A; the data indicates that £ ~h 2 j\. 

(iv) Intermediate: For h below the buckling threshold /ig, but above a crossover 
values which we denote by h T (see Figure 3A), a scaling of £ ~ h^ 5 is observed. 
Unlike the buckling threshold hg, the crossover thickness hj strongly depends on 
A. We note that an h^ 5 energy scaling has also been observed experimentally in 
[11] for the case of a single-stepped metric, even though the theoretical prediction 
in this work was slightly different (we will return to this point in the Discussion). 

Figure 3B displays a log-log plot of a mean longitudinal curvature k as function 
of the thickness h for the same parameters. Here too we obtain four distinct regimes: 

(i) Planar configuration: For h > h B the equilibrium configuration is flat and 
therefore k = 0. 

(ii) Near threshold: Just below the buckling threshold k is almost independent of 
the transition width A (Figure 3B, inset). 

(iii) Isometric configuration: As h — > the curvature tends to a constant, as expected 
by the prediction that the equilibrium configuration converges to a bending 
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Figure 3. (A) Total elastic energy £ (expressed in units of the stretching modulus 
Yh/8) versus h for A = 0.1 and three values of the transition width A (see legend). 
The solid curves correspond to the theoretical predictions E ~ h 2 and E ~ h 4 ^ 5 . (B) 
Mean longitudinal curvature k for the same parameters; the solid line corresponds to 
k ~ h~ 3/5 . The inset shows that there is almost no dependence of k on A close to the 
buckling threshold. 

minimizer among all isometries [28, 19]. The limiting curvature strongly depends 
on the transition width, k ~ 1/A. 

(iv) Intermediate: In the intermediate regime, h T «: h «: h B , we obtain a scaling of 

Figure 4A displays the function E(v) for A = 0.1, A = 1 /40, and several values of h. 
For small enough h, the value of E far from the transition zone is close to the reference 
values, 1 + A/2. The transition between the two bulk values is over a distance which 
we denote by Sg, and which varies with h. As h — > the transition width 5 E converges 
to the transition width, A, of the reference metric. 

In Figure 4B we plot an estimated value of the transition width 5e versus h for the 
same parameters as in Figure 3. The solid line corresponds to the theoretical prediction 
6 E ~ h^ 5 . In Figure 4C we plot an estimated value of the crossover thickness h T versus 
the transition width A. The solid line corresponds to theoretical prediction hj ~ h 5 ^. 

3.4. The isometric limit 

As h — > the equilibrium configurations converge to an isometric immersion of a 
that minimizes the bending energy [28, 19]. Thus, the first fundamental form of the 
limiting configuration is given by the reference values, 

E(p) = G(v) = f 2 (v). 
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Figure 4. (A) The metric component E versus v for A = 0.1, A = 1/40, and several 
values of h (see legend). (B) Log-log plot of an estimated value of 5e versus h for 
A = 0.1 and three different values of A (see legend); the solid line corresponds to the 
theoretical prediction 5e ~ h 4 ^ 5 . (C) An estimated value of the crossover thickness hj 
versus the transition width A; the solid line correspond to the theoretical prediction 
h T ~ A 5 ' 4 . 



Substituting E and G into (10), the second fundamental form is given by the one- 
parameter family of solutions: 



= V c / 2 -/' 



L = Jc f 2 - f' 2 and N = 



f' 2 -ff" 



and the corresponding bending energy is 



(f' 2 -ff") 2 \dv 

cf 2_ f ,2 j f 2- 
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The constant of integration c is determined by a minimization of £g. Differentiating 
£ B with resect to c we obtain the implicit equation: 

X w ( f'2 — ff") 2 

We have thus obtained a closed solution for the h — > limiting equilibrium 
configuration and the corresponding elastic energy. 

The solution depends on the transition width A through its dependence on the 
swelling profile /. Since L 2 > 0, it must be the case that for every v, 

cf\v)>r\v), 

which yield a lower bound for c. Substituting (2) and setting v = 0, it follows that 

A 2 



For small enough A and v in bulk regions we have / « 1 + A/2 and /' « 0, in which 
case 

£ ~ Vc(l+A/2)> A 

which indicates that the longitudinal curvature should grow (at least) proportionally 
to A" 1 ; the data in Figure 3B shows that this is in fact the correct scaling. 



3.5. The intermediate asymptotic regime 

We now turn to analyze the energy scaling in the intermediate regime, where h is 
below the buckling threshold h B/ but not too small (in a sense to be made precise). 

As seen in numerical computations, in this intermediate regime the metric a 
is very close to the reference metric a throughout most of the domain, except in a 
transition layer, which is larger than A. Moreover, the metric components E and G 
have different transition layers (data not shown), whose lengths we denote by 5 E and 
be, respectively; our computations show that <5g <^ b~E- Note that the composite sheet 
is mostly frustrated along the stripes, which means that the metric component G can 
more easily than to E approach its reference value. We will henceforth ignore the role 
of G in the energy minimization; a more detailed analysis (not shown) taking G into 
account reveals that this omission does not affect the results. 

We perform a scaling analysis to predict the dependence of the energy and the 
curvature on the various parameters. We start with the stretching energy £ s . The bulk 
values of the diagonal entries of the reference metric are 1 ± A/2, and the transition 
between these values is over a distance A. For A <sc <5 E , the metric discrepancy in E 
is of order A over a transition layer of order 5 E/ which yields a stretching energy of 
order 

£ S ~ A 2 6 E . (12) 
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We now turn to estimate the bending energy £ B . Consider (10): since the 
derivatives of E are small everywhere outside the transition layer, we have 

L 2 « c N ~ 0, for 17 not in transition zone, 

where we neglected A compared to 1. Inside the transition layer of size 5 E , E is of 
order one, whereas E' ~ A/5 E and E" ~ A/6 2 . Thus, inside the transition layer, 

, A 2 , A 2 

L ~ c r- and N ~ -—, (13) 

b\ L 2 5 E V ' 

and after straightforward algebraic manipulations the total bending energy (from 
both bulk and transition layer) satisfies the following asymptotic relation, 

E B ~h 2 wc + h 2 5 E (c-^] + h 2 5 E A " 



b\) ' ""-(c- A 2 /6|)6| 



, A 2 , , A 2 

Ifw^r + h z (w + <5 E )A + h L b E 



where A = c - A 2 /5 2 . We then further neglect 5 E compared to w in the middle term. 

The parameter A, which plays the same role as the integration constant c, does not 
affect the stretching energy, and therefore its value is determined by a minimization 
of the bending energy only. We get A ~ Aar 1/2 S~ 3/2 , which upon substitution back in 
the bending energy yields 

h 2 Aw(^ 

AW-- (14) 

We next turn to determine the value of the transition width 5 E . Unlike A, 5 E affects 
both the stretching and the bending, and is therefore determined by a minimization 
of the total energy, 

£ ~ A 2 6 ^ 




It is unclear a priori whether A dominates yfb^Jw or vice versa in the bracketed 
expression; both are small parameters. A straightforward asymptotic expansion 
reveals there are indeed two possible scalings for b E , 

w V3 h 2/3 if h/w «: A 3 (Case 1) 

A" 2/ V /5 /z 4/5 if h/w » A 3 (Case 2). 

Case 1 corresponds to a regime in which the bending is dominated by L 2 in the bulk, 
whereas Case 2 corresponds to a regime where the bending is determined by a balance 
between L 2 in the bulk and N 2 in the transition layer. 

These two asymptotic regimes yield the following scalings for the energy, 

A 2 w 1/3 h 2/3 if h/w « A 3 (Case 1) 
A 8/5 w V5 h 4/5 if h/w » A 3 (Case 2), 



Metrically discontinuous non-Euclidean plates 



13 



and the following scalings for the longitudinal curvature, 

A 1/2 wT 1/2 /T 1/2 if h/w « A 3 (Case 1) 
A 4/5 aT 2/5 r 3/5 if h/w » A 3 (Case 2). 

A comparison to our computational results indicates that we obtained scaling 
laws that match Case 2. Indeed, A 3 ~ 10" 3 , which means that for Case 1 to occur 
we need h/w <sc 10" 3 , which is already close to the isometric regime. It takes simple 
algebra to see that three distinct regimes are obtained only if A/w <sc A 2 , in which 
case, 

' A 2 wA~ 2 h 2 if h «: w(A/w) 3/2 (Isometric limit) 
A 2 w 1/3 h 2/3 if w(A/w) 3/2 « h « wA 3 (Case 1) 
A 8/5 w 1/5 h 4/5 if wA 3 « h « h B (Case 2). 

If, however A/w <£. A 2 , then we obtain only two distinct regimes, 

A 2 wA~ 2 h 2 if h «: h T (Isometric limit) 

A 8/ V /5 fr 4/5 if hj « h/w « h B (Case 2), 



where 



/A\ 5/4 
h T = A 1 ' 2 ™ I £j 



In Figure 4C we plot the dependence of for on A and find that the predicted scaling is 
indeed correct. 



4. Discussion 



This paper presents a novel analysis of thin sheets of pre-stressed materials, in which 
the reference metric has a uniaxially periodic structure. Such a geometry is common 
in nature (e.g. in monolayers of fibrous tissue) as well as in man-made composite 
materials (e.g., [1]). As explained in Subsection 3.1, such an incompatible metric 
favors a combination of longitudinal bending and transversal undulations. It should 
be noted, however, that the occurrence of other strain relieving mechanisms, such as 
wrinkles has not been ruled out. 

Another characteristic of the reference metric is the abrupt transition between 
two locally Euclidean metrics. We showed that in such setting there is no simple 
distinction between a "thick plate" and a "thin plate". When varying the thickness 
of the sheet, the equilibrium conformation switches between several regimes, each 
satisfying its own scaling laws, thus reflecting different balances between energetic 
contributions. 

Our results are closely related to the theory and experiments reported in [11]. 
Rather than considering an array of alternating stripes, they considered a system 
which in our context corresponds to a single period. In addition, their experimental 
setting corresponds to a regime of parameters in which A «a), i.e., a very abrupt 
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metric transition, or equivalently, wide stripes. Their theoretic predictions were 
scaling relations of £ ~ h 2/3 and k ~ h~ 1/2 , which coincide with the regime we labeled 
Case 1. Yet, their numerical calculations clearly exhibit an energy scaling of E ~ h 4/5 , 
(Figure 2G in [11]) which rather corresponds to our Case 2 regime. Thus, it appears 
as if their parameters correspond to the latter regime, which was only discovered in 
the present work. It should be noted, however, that £ ~ h 2/3 is indeed the generic 
scaling for sufficiently small A and h. 

The revelation that distinct scaling regimes may emerge with decreasing 
thickness is the most important contribution of this paper. We conjecture that it 
is a general feature in non-Euclidean plates in which the metric has a multi-scaled 
structure. At finite thickness h, bending dominates all features below a certain length 
scale, so that metric discrepancies below this scale will remain unnoticed. Effectively, 
this is equivalent to the metric being smoothed at a cutoff scale determined by the 
thickness. As h is further reduced, this cutoff scale decreases, and new features of the 
metric may thus manifest. 
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